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INTRODUCTION 


In  support  o£  an  Army- wide  modernization  program,  the  Army 
Armament  Research  and  Development  Command,  Dover,  NJ,  is  engaged  in 
the  development  of  safety  criteria  as  a project  entitled,  ’^Safety 
Engineering  in  Support  of  Ammunition  Plants.'* 

Prior  to  investigating  the  explosive  fragmentation  of  reinforced 
concrete  structures  (one  aspect  of  the  safety  engineering  program) , 
a feasibility  study  of  sub-scale  dynamic  modeling  was  conducted  to 
determine: 

1.  The  applicability  of  sub-scale  modeling  to  random  mass 
fragments  of  reinforced  concrete  subjected  to  excessive  blast  loads. 


2.  The  physical  relationships  between  the  model  and  the  pro- 
totype in  constitutive  properties,  panel  geometry,  and  boundary  con- 
ditions as  well  as  the  resulting  fragment  weight/shape,  density,  and 
distribution. 


PROBLEM  DEFINITION  AND  SOLUTION 


In  munitions  manufacturing  facilities,  reinforced  concrete 
plays  an  essential  role  as  a building  material  for  both  dividing 
walls  (between  quantities  of  explosive)  and  exterior  walls  and 
roofs. 

Mechanics  of  Fragmentation  (ref.  1) 

When  a reinforced  concrete  element  is  substantially  overloaded 
by  a blast  wave,  the  element  fails  and  concrete  fragments  (post- 
failure) are  displaced  at  high  velocities  (fig.  1). 

Failure  of  an  element  is  characterized  by  the  dispersal  of  con- 
crete fragments  formed  by  the  cracking  and  displacement  of  the  con- 
crete between  the  donor  and  receiver  layers  of  the  reinforcement. 

As  an  element  deflects  and  the  concrete  begins  to  crush,  the  compres- 
sion stresses  normally  resisted  by  the  concrete  are  transferred  to 
the  reinforcement.  With  increased  deflections,  these  compression 
forces  tend  to  buckle  the  reinforcement  outward,  initiating  rapid 
disintegration  of  the  element. 
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The  velocity  o£  individual  fragments  varies  and  depends  on: 

(1)  the  magnitude  of  the  excess  (or  blast)  impulse  minus  the  flexural 
impulse  capacity  of  the  element  (area  under  the  resistance-time 
curve),  (2)  the  mass  of  the  fragment,  (3)  the  location  of  the  frag- 
ment prior  to  collapse,  (4)  the  interaction  between  fragments  during 
their  flight,  and  (5)  the  strength  and  time  history  of  the  compres- 
sive stress  wave  transmitted  through  the  wall  as  the  blast  wave  is 
reflected.  Although  the  velocities  of  individual  fragments  differ, 
the  average  translational  velocity,  v^  (avg),  of  the  debris  after 
complete  failure  can  be  approximated  from  the  excess  impulse,  i^,  and 
the  unit  mass,  m,  of  the  barrier;  i.e.,  the  momentm  of  the  wall 
after  collapse  is  numerically  equal  to  the  excess  impulse,  i^  = mv^ 
(avg) . 


H 

where 

i = applied  unit  blast  impulse 

pH  = reinforcement  ratio  in  the  horizontal  direction 

d^  = distance  between  the  centroids  of  the  compression 
and  tension  reinforcement. 

f^^  = dynamic  design  stress  for  the  reinforcement 

H = span  height 

V£  = maximum  velocity  of  the  post-failure  fragments 
= impulse  coefficient  for  ultimate  deflection 
= post-failure  fragment  coefficient 

The  reflected  shock  wave  is  also  transmitted  through  the  wall 
as  a compressive  wave,  and  reflects  from  the  rear  surface  as  a ten- 
sion wave  which  partially  cancels  the  incoming  compressive  wave.  If 
the  shock  strength  is  great  enough  and  its  decay  rate  rapid  enough, 
a net  tensile  stress  will  occur  at  some  distance  from  the  rear  sur- 
face. The  concrete  (which  is  weak  in  tension)  will  fail  and  spall 
off.  The  thickness  of  the  spall  and  its  velocity  can  be  estimated 
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from  knowledge  of  the  reflected  wave  properties  and  the  constitutive 
properties  of  the  concrete. 

Dynamic  Modeling  of  Reinforced  Concrete  Elements 

The  following  exercise  defines  the  physical  similarity  between 
reinforced  concrete  elements  in  the  model  and  the  prototype  subjected 
to  blast  overloading  as  described  in  reference  2.  The  constitutive, 
geometric,  and  mechanical  parameters  dictated  by  the  scope  of  the 
problem  are  given  in  table  1. 


PI  TERMS  FROM  DIMENSIONAL  HOMOGENEITY 

By  inspection  (already  dimensionless) 

= K 

Technique  from  reference  2:  Matrix  Solution 

pOLO.j.od  34^  a5^  a0y  ay^  a3^  ag^  ^IOt  ^^2 

csccsscs  A s 

E ^13w  ai4^ai5 

c F 

Using  matrix  algebra  (identity  submatrix) 

Form  a matrix: 
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R p. 

Ps 

c 

c 

T 

c 

u 

s 

T 

s 

L 

c 

L 

s 

E E Wt, 

s c F 

V 

a 

a a 

1 2 3 

a 

4 

a 

5 , 

a 

a 

7 
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ft 
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ft 
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12  13  14 
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F 
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+1 
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0 

0 2 
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0 
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Matrix  algebra  steps; 

Step  1 - Divide  row  3 by  row  2. 

Step  2 - Subtract  new  row  3 from  row  1,  and  write  new  row  1. 

Step  3 - Subtract  row  1 from  row  2. 

Step  4 - Multiply  row  3 by  row  4 and  add  to  row  2.  Write  new 

row  2 . 

Therefore,  identity  submatrix  in  the  following  matrix; 


a 

1 a2 

as  a4 

as 

as 

a ' 

7 a0 

ag 

a 1 0 

a 

1 1 a 

12  a 1 ; 

3 a 

14  a 1 5 

F 
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0 

0 

0 

0 

0 

0 

0 

0 

0 
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Then, 

by  in; 

spection 

ai 

= 

-as 

-as  -ay 

-as 

-ai  2 

-a 

1 3 -a 

1 4 

az 

= 

335 

+ 3as  + 
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+ 3a 

8 - 

ag  - 

a 1 0 

-ail 

+ 

3ai2 

+ 3a 

1 3 

-3ai  5 
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Substituting  in  the  original  equation 

f0l0t0=  w'^S-ae-ay-as-ai  2-ai  3-ai  It  r ^^5'*'^^6'^^^7‘^^^8"^9"^10"^11 
+3ai 2+3ai 3-3ai 5 

X P ^IOt  ail  ai2 

a a c s s c s 
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Collecting  terms 
be  expressed  as: 

of  like  exponent,  we  have 

the  model 

law.  This  can 

, \ 

W 

^ P 

r3c  R^T  R^U  R^T 

L L 

L.  R^E  r3e 

V 

cess 

c s 

A s c 

r3 

W W , W , W 

, R , R , 

R , W , W 

jL 

K 

J 


The  law  can  also  be  stated  in  the  form  of  table  2. 


SCALE  FACTORS 


Various  relations  between  scale  factors  can  be  determined  from 
the  dimensionless  terms. 


TTi  X = A 


= 1 (set  equal  to  1 for  ease  in  fabricating  model) 


"^7 


3 

A ^ R C C S 5 


= x„  = X, 


'^9>  '^10  ^ ^ ^ 

c s A 


•^11 

TT  I 2 
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These  interrelations  are  all  satisified  by  a replica  model ^ 
with  geometry  identical  to  the  prototype.  Although  wave  transmission 
effects  are  not  shown  in  scale,  they  will  scale  properly  because  all 
constitutive  properties  and  densities  are  unchanged. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  preceding  similarity  exercise  has  demonstrated  that 

1.  A replica  model  accurately  reflects  the  performance  of  a 
prototype  system. 

2.  Failure  modes  of  a replica  model  are  representative  for 
geometrically  similar  panels. 

3.  The  model  also  exhibits  similar  relationships  for  frag- 
ment mass  and  velocity,  distribution,  and  debris  density. 

For  economic  reasons  and  ease  in  testing,  it  is  recommended 
that  subscale  replica  models  be  used  with  a scale  factor  as  small  as 
possible  without  sacrificing  the  quality  of  the  model. 
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^A  replica  model  is  a physical  model  of  a prototype  which  is  geo- 
metrically similar  in  all  respects  and  employs  the  same  materials 
in  the  same  locations  as  the  prototype  (ref.  2). 
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3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 


Table  1.  List  of  parameters  for  explosive 
fragmentation  of  reinforced 
concrete  elements 


Parameter 

w 

R 

Pc 

Ps 


Dimension 

FL 

L 

FT^/L** 

FT^/L** 

F/L^ 


U, 


F/L^ 

F/L^ 

F/L^ 

L 


W. 


L 

F/L^ 

F/L^ 

FL 


V 


Description 

Energy  in  blast  source 

Radial  distance  of  blast 

Density  concrete 

Density  steel 

Compressive  strength 
(concrete) 

Tensile  strength  (concrete) 

Ultimate  strength  (steel) 

Tensile  strength  (steel ) 

Characteristic  geom . of 
element 

Characteristic  geom . of 
rebar 

Aggregate  size 

Elastic  moduli  (steel) 

Elastic  moduli  (concrete) 

Energy  imparted  into 
fragments 

Vol  of  fragments 

Distribution  function  of 
fragments 

% energy  absorbed  by 
element 
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Table  2.  PI  terms  - explosive  fragmentation 


Tfl 
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TT2 

= 
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= 
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TT^ 

= 
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= 
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Try 
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> 

ITS 

= 

Lc/R 

TTg 

= 

Ls/R 

> 

GEOMETRIC  SIMILARITY 

■^10 

= 

LA/R 

> 

W|,/W 

SIMILAR  ENERGY  TRANSFERS 

1712 

= 

V/r3 

SIMILAR  MASS  TRANSFERS 

■^13 

= 

DISTRIBUTION  FUNCTION 

iti4 

= 

K 

^ ENERGY  ABSORBED 
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REINFORCED  CONCRETE 

donor  protective  element 

(Explosive  Charge) 


LETHAL  ZONE 


ENERGY  AND  IMPULSE  BALANCES 


Blast  Energy  In  = Strain  Energy  + Sum  of  Fragment  Kinetic  Energies 
= Sum  of  Fragment  Kinetic  Energies 

(]  - K) (where  K = % Blast  Energy  Absorbed) 

Applied  Impulse  = Flexural  Impulse  Capacity  + Excess  Impulse 


ENERGY  IN 


BLAST 

LOADING 


FRAGMENT  ENERGY 


M3V3 


MjVi 


POST  DETONATION 


ENERGY  ABSORBED 


Figure  1.  Explosive  fragmentation. 
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